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We show that the lattice calculation of K — » nn and e'/e ampli- 
tudes for (8,1) and (27,1) operators to 0(p 4 ) in chiral perturbation 
H-j , theory is feasible when one uses K — > hit computations at the two 

\Q ' unphysical kinematics allowed by the Maiani- Testa theorem, along 

01 , with the usual (computable) two- and three-point functions, namely 

K —* 0, K — > 7r (with momentum) and K-K. Explicit expressions 
, for the finite logarithms emerging from our 0{p i ) analysis to the 

' above amplitudes are also given. 

o \ 

cr> ' 1 Introduction and Motivation 

o : „ 

CN ■ Recent lattice QCD calculations done by the CP-PACS @ and RBC g collabo- 

rations using domain wall fermions have made significant progress in explaining 
r"| ■ the AI — 1/2 rule in the decay K — ► 7T7t, though their results for the direct CP 

Oh! violation parameterized by Re(e'/e) differs rather drastically from experiment. 

Recall that measurements at CERN j|] and Fermilab Q have yielded an ex- 
\ perimental grand average of Re(e'/e) = (17.2 ± 1.8) x 10 -4 . On the theoretical 

' side, both lattice collaborations find a value of Re(e'/e) ~ —5 x 10 -4 , a negative 

value, though both groups have made rather severe (uncontrolled) approxima- 
tions. Given that large cancellations occur between contributions of the strong 
and the electroweak penguins towards e'/e (cancellations that are not relevant 
jj^ ■ in the calculation of the CP-conserving K — » tttt amplitudes), and given the 

serious approximations, the disagreement with experiment for e'/e should not 
be totally unexpected ||, 

One of these uncontrolled approximations was the use of the quenched ap- 
proximation, where the fermion determinant in the path integral is set to 1 in 
order to make the problem tractable on current computers. Another was the 
use of leading order chiral perturbation theory to relate unphysical K — > n and 
K — > |0) amplitudes to the physical K — > itit amplitudes. This method was first 
proposed by Bernard, et al. 0. Because of the difficulty of extracting multi- 
hadron decay amplitudes from the lattice, as expressed by the Maiani- Testa 
theorem it is much easier to compute the two- and three-point functions 
(i.e. K — > |0) and K — » ir, respectively) and use chiral perturbation theory 
(ChPT) to extrapolate to the physical matrix elements. 
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It is likely, however, that next-to-leading order ChPT will introduce signifi- 
cant corrections (~ 30% or more) to the leading order amplitudes. Furthermore, 
since final state (strong) phases cannot arise at tree-level in the chiral ampli- 
tudes, chiral-loop corrections are essential to enable us to use the measured 
phases for the 1 = and 2 final states as additional testing ground of the 
calculational apparatus. Unfortunately, at higher orders in ChPT the number 
of free parameters that enter the theory (and must be determined from first 
principle methods like the lattice) proliferates rapidly. It has been shown by 
Cirigliano and Golowich 0] that the dominant electroweak penguin contribu- 
tions [ (8l,8r)'s ] to K — + tttt can be recovered at next-to-leading order (NLO) 
from K — > 7r amplitudes using momentum insertion. Bijnens et al. ficfl showed 
how to obtain most of the low energy constants (LEC's) relevant for the case of 
the (8l, 1r)'s and (27^, 1r)'s using off-shell K — > tt Green's functions; not all 
LEC's could be determined, however. 

On the lattice, though, not only K — > |0) and K — ► tt with momentum 
insertion are calculable, but so is K — > tttt at the two values of unphysical kine- 
matics for which the Maiani- Testa theorem can be bypassed. To recapitulate, 
despite Maiani- Testa restrictions, direct calculation of K — > tttt on the lattice is 
accessible at (i) m^ ttICC = rn}^ ttlcc where the weak operator inserts energy pi 
and (ii) m^ tticc = 2m^ attico , i.e. at threshold @. We will refer to these two 
special locations as unphysical kinematics point 1 (UK1) and point 2 (UK2), 
respectively. In this work we therefore focus on using K — ► |0), K — > tt with 
momentum insertion and K-K, along with information from K — > tttt at these 
two unphysical values of the kinematics which are accessible to the lattice Jl3|. 
Thereby, we are able to show that all the relevant 0{p A ) LEC's can be recov- 
ered for K — > 7T7T in the physical (8l, 1r) and (27l, 1r) cases. Expressions for 
0(p 4 ) finite logarithimic contributions to all the processes that may be needed 
for fitting the lattice data are then given. 

The content of the paper is as follows. Section 2 very briefly recapitulates the 
formalism of effective four-fermion operators in a Standard Model calculation. 
Section 3 reviews ChPT and the realization of the effective four-quark operators 
in terms of ChPT operators for weak processes. Section 4 presents the results 
of this paper, showing how to obtain the low energy constants necessary for 
physical K — ► tttt amplitudes at 0{p A ) from quantities which can, in principle, 
be computed on the lattice. Section 5 presents the conclusion. Finally, the finite 
logarithmic contributions to the relevant amplitudes are presented in a set of 
Appendixes. 



2 Effective Four Quark Operators 

In the Standard Model, the nonleptonic interactions can be expressed in terms of 
an effective AS = 1 hamiltonian using the operator product expansion [l5| , 

<vrvr|WA5=i|^) = ^V ud V: s Y^^){^\Q t \K)„ (1) 
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where Ci([i) are the Wilson coefficients containing the short distance perturba- 
tive physics, and the matrix elements {k^Q^K) ^ must be calculated nonpcr- 
turbatively. The four quark operators are 



Q 


1 = Sa7p( 1 


- J 5 )u a u b ^(l - 7 5 )4, 


(2) 


Q 


2 = S a 7p(l 


- 7 5 )u 6 u & 7 fl (l - 7 5 )da, 


(3) 


Q 3 


= Sq7 m (1 - 


7 5 )da^g fc 7 P (l-7 5 )9 6 , 

9 


(4) 


Qa 


= s q 7m( 1 - 


7 5 HE^(l-7 5 )9a, 

9 


(5) 


Q 5 


= s q 7m( 1 ~ 


7 5 K^9 6 7 P (l+7 5 )%, 
9 


(6) 


Q 6 


= ^7^(1 - 


7 5 H^ 6 7 AI (l+7 5 )<Za, 
9 


(7) 


Qi = 




7 5 K^e 9 g fc7 ^(l + 7 5 )g 6 , 
9 


(8) 


Qs = 




7 5 H^e g g 6 7' t (l + 7 5 )9a, 
9 


(9) 


Q9 = 


3 - n 


7 5 K^e 9 g fc 7^(l-7 5 )gfc, 
9 


(10) 


Q10 = 




-7 5 )4^e^ 6 7' l (l-7 5 )g a . 


(11) 



In the effective theory Q\ and Qi are the current-current weak operators, 
Q3 — Q& are the operators arising from QCD penguin diagrams, while Q-j — Q10 
are the operators arising from electroweak penguin diagrams. 

3 Chiral Perturbation Theory 

Chiral perturbation theory (ChPT) is an effective quantum field theory where 
the quark and gluon degrees of freedom have been integrated out, and is ex- 
pressed only in terms of the lowest mass pseudoscalar mesons ]l6| . It is a 
perturbative expansion about small quark masses and small momentum of the 
low mass pseudoscalars. The effective Lagrangian is made up of complicated 
nonlinear functions of the pseudoscalar fields, and is nonrenormalizable, making 
it necessary to introduce arbitrary constants at each order in perturbation the- 
ory. In such an expansion, operators of higher order in the momentum (terms 
with increasing numbers of derivatives) or mass appear at higher order in the 
perturbative expansion. The most general set of operators at a given order can 
be constructed out of the unitary chiral matrix field S, given by 
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E = cxp 



~2icj) a \ a 



(12) 



where A a are proportional to the Gell-Mann matrices with tr(A a Afc) = 6 a b, <fi a 
are the real pseudoscalar-meson fields, and / is the meson decay constant in the 
chiral limit, with equal to 130 MeV in our convention. 

At leading order [0(p 2 )] in ChPT, the strong Lagrangian is given by 



'-st 



tr^Ed^E] 



f 2 B 



(13) 



where \ = ( m u, m d , m s ) diag and 
Bo 



T + 



The leading order weak chiral Lagrangian is given by [p 



= a 8 8tr[A 6 EQEt] + aitr[A 6 S |t E^St] + a 2 2B tr[A 6 ( x t E + E+ X )] 
+a 27 q i (m^) k i {m^)) + H.c, (14) 

where is symmetric in i,j and traceless on any pair of upper and lower 
indices with nonzero elements = lj ^22 = 1/2 and ijjf = —3/2. Also, Q is 
the quark charge matrix, Q — 1/3(2, —1, — l)diag and (Ae)ij — ^3^2- 

The terms in the weak Lagrangian can be classified according to their chi- 
ral transformation properties under SU(3)l x SU(3)_r. The first term in (14) 
transforms as 8l x 8_r under chiral rotations and corresponds to the electroweak 
penguin operators Q7 and Q$- The next two terms in (14) transform as 8l x 1r, 
while the last transforms as 27^ x 1r under chiral rotations. All ten of the four 
quark operators of the effective weak Lagrangian have a realization in the chiral 
Lagrangian differing only in their transformation properties and the values of 
the low energy constants which contain the non-perturbative dynamics of the 
theory. 

For the transition of interest, K — > nn, the operators can induce a change 
in isospin of ^ or | depending on the final isospin state of the pions. We can 
then classify the isospin components of the four quark operators according to 
their transformation properties & ^ : 



l /2 > 




?9 /2 ,Q^ 2 :8 L x 1«©27 L x 1 R 


3/2 


Q2 > £ 


f 9 /2 ,Q% 2 : 27 L xl R ; 






Q\ /2 ,Ql /2 ,Ql /2 -SlxIu; 




Qi 1 
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Note that Q3 — Qe are pure isospin \ operators. This paper deals only 
with the 27 l xIr and 8l X 1r operators. For the treatment of the 81 x 8jj 
operators to 0(p 2 ) NLO see Ref. ||. At NLO the strong Lagrangian involves 
12 additional operators with undetermined coefficients. These were introduced 
by Gasser and Leutwyler in jl8| . The complete basis of counterterm operators 
for the weak interactions with AS" =1,2 was treated by Kambor, Missimer and 
Wyler in Jl9| and p0| . A minimal set of counterterm operators contributing to 
K —* it and K irir for the (8l, 1j?) and (27l, 1r) cases is given by |l7| , with 
the effective Lagrangian 

c w = Yl e ^ (8,1) + J2 d M 27,1 \ (15) 



(8,1) 
0(8,1) 
0(8,1) 



(8,1) 
4 

(8,1) 
5 

(8,1) 
10 

/n(8,l) 
c/ 12 



(8,1) 
13 

(8,1) 
15 

(8,1) 
35 

(8,1) 
39 



tr[A 6 S 2 ], 

tr[A 6 S]tr[S], 

tr[A 6 P 2 ], 

tr[A 6 P]tr[P], 

tr[A 6 [5,P]], 

tr[A 6 {S,L 2 }], 

ti[X e L^SL% 

tr[AeL„]tr[{L",S}] 1 

tr[A 6 S][L 2 ], 

tr[A 6 [P,L 2 ]], 

tT\X 6 {L^d v W^}], 

tx[X e W^W' 1 % 



0(27,1) 
0(27,1) 



(27,1) _ 
4 — 
(27,1) 
5 

(27,1) 



0(27,1) 



o 
a 



(27,1) 
20 

(27,1) 
24 



4(L^(L^tT[S], 

t%(L^(d„wn l v 



(16) 



with S = 2P (x f S + St x ), P = 2P (x f S - E+x). L n = iE^E , and W**" = 
2(d„L v + 

This list is identical to that of Bijnens et al. jlO| except for the inclusion 
of 035' 3g and 24 which contain surface terms, and so cannot be absorbed 
into the other constants for processes which do not conserve 4-momentum at 
the weak vertex. Since we must use 4-momentum insertion in a number of our 
amplitudes, these counterterms must be considered, and they are left explicit 
even in the physical amplitudes. There are additional operators containing 
surface terms, but it was checked that these counterterms can be absorbed into 
linear combinations of the above minimal set for all amplitudes considered in 
this paper. 

The AS = 2 operators are components of the same irreducible tensor pl| 
under SU(3)l x SU(3).r, and so the d% are the same for both the AS = 1 
and AS = 2 cases. The operators governing AS* = 2 transitions are obtained 
from the above (27^, 1.r)'s, only with — i 2 | = ljife/ = otherwise. This is 
important since some of our information comes from the K° — > K° amplitude. 
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The divergences associated with the counterterms have been obtained in jl9| 
and JlQ[. The subtraction procedure can be defined as 



1 



16tt 2 / 2 



2(a 1 e i + a 2 e' i ), (17) 



16tt 2 / 2 



:{lE - 1 - ln47r) 



(18) 



with the divergent pieces, £,,£^,7, given in Table 1. 

It is also necessary for the method of this paper to consider the 0(p A ) strong 
Lagrangian, which was first given by Gasser and Leutwyler, = Y] LjO^ . 



Table 1: The divergences in the weak 0(p A ) counterterms, e^'s and di's, for the 
(8,l)'s and (27,l)'s, respectively. 





£j 


< 


di 


li 


1 


1/4 


5/6 


1 


-1/6 


2 


-13/18 


11/18 


2 





3 


5/12 





4 


3 


4 


-5/36 





5 


1 


5 





5/12 


6 


-3/2 


10 


19/24 


3/4 


7 


1 


11 


3/4 





20 


1/2 


12 


1/8 





24 


1/8 


13 


-7/8 


1/2 






15 


23/24 


-3/4 






35 


-3/8 









39 


-3/16 










The strong 0(p 4 ) operators relevant for this calculation are the following Esl : 



(st) 


-tr[L 2 ] 2 , 


(st) 


= tr^i^tr^L^ 


(st) 


= tr[L 2 L 2 ], 


(st) 


= tr[L 2 ]tr[S], 


Q (st) 


= tr[L 2 S], 


n (st) 
^6 


= tr[5] 2 , 


n {st) 
^8 


= \tv[S 2 -P 2 ]. 



The Gasser-Leuytwyler counterterms also contribute to the cancellation of 
divergencies in the expressions relevant to this paper. The subtraction is defined 
similarly to that of the weak counterterms, 
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U = I.: 



16tt 2 



1 1/ 

^4 + 2 (7£ 



1 -ln47r) 



(20) 



with the divergent parts of the counterterm coefficients given in Table 2 |p 



Table 2: The divergences in the strong 0(p 4 ) counterterms, r^Jl8[. 



i 


r, 


1 


3/32 


2 


3/16 


3 





4 


1/8 


5 


3/8 


6 


11/144 


8 


5/48 



4 K — > 7T7T amplitudes at 0(j? 4 ) 

As mentioned before, in this work we will include both K — > 7r with momentum 
insertion and A' — * 7T7t at the two unphysical kinematics. The complete list of 
necessary ingredients consists of: the two point functions K° — > |0), the three 
point functions K° — > A , and K — > n, all with m s ^ m<i = m„, and the four 
point functions K — > 7T7t at the two values of unphysical kinematics, m^- = m w 
(requiring energy insertion), and mj( = 2m^. These two threshold values of the 
kinematics bypass the Miani- Testa theorem, which states that multihadron final 
states are not accessible on the lattice at any other kinematics aside from the 
threshold Q. At these kinematics, the strong phases are 0, and the effects of 
final state interactions vanish. However, these amplitudes at unphysical kine- 
matics do contain information on the 0(p 4 ) low energy constants, and when 
combined with information from the other two- and three-point functions men- 
tioned above, all of the 0(p A ) low energy constants necessary for K — > tttt can be 
obtained. The phases of the amplitude are introduced in ChPT via the one-loop 
unitarity corrections of the 0(p 2 ) operators. 

Because the K — > 7r amplitudes do not conserve four-momentum for m s =/= 
rrid, it is necessary to allow the weak operator to transfer a four- momentum 
q = pk — P-k i as in || . This is also necessary for the case of K — > tttt, wlk = rn^ 
@. At 0(p 4 ), this requires the inclusion of (potentially many) surface terms 
in our minimal counterterm operator basis. The number of additional such 
terms appearing in linearly independent combinations was discovered to be small 
(four), and thus the method was not invalidated. Also, this method requires 
the computation of K — + tttt matrix elements at unphysical kinematics because 



7 



there are LEC's which appear in K — > irir but do not appear in K —* 7r at all. 
These are c?5,ei3 and ei5. 

The diagrams which must be evaluated for the 0(p 4 ) corrections are shown 
in Figure 1. The diagrams to be evaluated for K — > |0) are Al and A2. The 
diagrams to be evaluated for K — > 7r and K ^ K are B1-B3. C1-C6 and D1-D6 
must be evaluated for if — ► 7T7t. D1-D6 contain the tadpole vertex of the weak 
mass 0(p 2 ), (8,1) operator. Also, the renormalization of the external legs via 
the strong interaction must be taken into account. 

4.1 (27 L ,l fl ), A/ = 3/2 

The counterterms necessary to reconstruct the (3(p 4 )(27l, 1r), Ai = 3/2, K — > 
7T7T amplitudes can be obtained from if — ► if ; i4T + — > 7r + ,Ai = 3/2; and 
if — > 7T7T, Ai = 3/2 at either value of the unphysical kinematics. The expres- 
sion for if — > if is given by (all masses and decay constants are the bare ones) 



(K \O^\K°U = -^m 2 K -j^[4(d[ + d r 2 + d r 20 -4d r 24 -cf 4 -d r 7 )m% 
-2(44 + d r 7 )m 2 K ml + 4d£ro*] (21) 

Equation (21), as well as all the following amplitudes, include only the tree 
level 0(p 2 ) and 0(p 4 ) weak countcrtcrm contributions. For brevity, the logarith- 
mic terms, as well as the Gasser-Leutwyler L t counterterms have been omitted 
in the body of the paper, but are included in a set of Appendixes. It was verified 
that the divergences in the logarithmic terms cancel those of the counterterms, 
providing a strong check on the calculation. Note also that for the application 
of this method most of the Gasser-Leutwyler counterterms must be known, and 
that an improved determination of the relevant ones could be obtained from a 
lattice calculation of observables in the purely strong sector, e.g. most can be 
obtained from the pseudoscalar masses and decay constants. 

From the above K° — > if amplitude, one can extract the values of d[ and d 7 
from a fit to terms quadratic in the quark masses. The other relevant expressions 
for K -> 7T7T, Ai = 3/2 are 



(7r +| (27,i),(3/2) |x+)ct = _jL a27pK . pw + j-[2d r 2m 2 K ml + (d r 20 - d\ - 2d r 7 ) 
m 2 K p K ■ Ptt + (d 20 -d r 4 - d r 7 )m 2 p K ■ Pn ~ 8d 24 (p K ■ Pn) 2 ], (22) 



( 7r + 7r -| O (2r,l),(3/2)| jr >ct = _!^ m 2 + 1^! 



(d£o - 2dl + d r 5 - 3d?), (23) 



for K — > 7T7T, rriK — — m, and 
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Al A2 Bl B2 




B3 CI C2 C3 




C4 C5 C6 Dl 




D2 D3 D4 D5 




D6 

Figure 1: 0(p 4 ) corrections include tree-level diagrams with insertion of the 
0(p 4 ) weak vertices (crossed circles), tree-level diagrams with insertion of 0(p 4 ) 
strong vertices (lightly shaded circles), one-loop diagrams with insertions of the 
0(p 2 ) weak vertices (small filled circles) and the 0(p 2 ) strong vertices (big filled 
circles). Al and A2 are for K — » |0). B1-B3 arc for K -> tt and K —>~K. C1-C6 
and D1-D6 are for K — > nir. 
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(24) 

for K -> 7T7r,TO K(1 _ i00 p) = 2m 7r ( 1 _ ioop ). 

From Eq. (22) we get the additional combinations of counterterms d 2l d 2i , 
and d\ — d 20 . From either expression for K — > tttt at unphysical kinematics we 
can then obtain d\ — cfij. Along with the tree level LEC, «27, these five linear 
combinations [d^d^, d% — d£, dj— dJjoi ^24] are sufficient to determine -FT — > tttt at 
the physical kinematics as given in Eq. (25). Notice that there is considerable 
redundancy in determining these coefficients. For example, d± — d2o, d^ — d$, 
d 2 ... occur in several of Eqs. (21)-(24)|||. 



(7r+7r -|O( 2 T > i),(3/2) | ^0 )ct = -^ {m \ _ ml)^ loop) +j^(m%- ml) 
[(-dj + d r b - 4d r 7 )m 2 K + (44 + AdT 20 - 16d r 2i - 4d£ - 2d r 7 )ml] . (25) 

The logarithmic and Gasser-Leutwyler counterterm contributions to the ex- 
pressions in this subsection are given in Appendix C. Note, also, that for the 
cases of physical K ^ tttt amplitudes, Eqs. (25), (34), and (35), the pseudoscalar 
decay constants and masses are the physical (renormalized to one-loop order) 
ones. For all other amplitudes given in this paper except K — ► tttt at physical 
kinematics the formulas are in terms of the bare constants. The distinction 
between bare and renormalized constants is made only in tree-level amplitudes, 
since making this distinction in the 0(p 4 ) expressions introduces corrections at 
higher order [0(p 6 )] than is considered in this paper. 



4.2 (8 L ,l fl ) + (27 L ,l fl ), A/= 1/2 

The counterterms necessary to reconstruct the 0(p 4 ) [(8^, 1 r)+{27l, 1r)], AI = 

1/2 

1/2, K — > tttt amplitudes, relevant for operators such as Q 2 , can be obtained 
using di's obtained from the [(27, 1); AI = 3/2] case given above along with 
information from K° -> \0);K+ -> tt+,AI = 1/2; and K -> tttt, AI = 1/2 at 
both unphysical kinematics. For K° — > |0), we have 



(0|O (8 



'-(m 



K 



/ 



r \ 

e 5 )m 



+ (2el + e r 2 - 2e r 5 )m 2 K ml + e r 2 m\\ 



(26) 



{V\0^\K») ct = -^d\{ m l- m lf. (27) 

Given the previously obtained value of d\ from the AI = | case, we can 
obtain e 2 and e\ — e£ from K° — > |0). The other relevant expressions are 
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( 7r +|0(27,l),(l/2)|^+ )ct = -JL a27pK . ?7r _ A[ 6d r m ^ _ 2 (3dJ + d^m* 

+ (-4o + ^4 - 3cJ 6 + 2d r 7 )m 2 K p K ■ p v + (-d 2Q + dl + 3rfg + d 7 )m\p K ■ p n 

+ 8d r 24 (p K ■ Pit) 2 }, (28) 



(7r +| (8,i),(i/2) |^+ )ct = i_ aim . ^ _ i_ a2 „& _ ^[ 2 ( e ; + e r - e g) TO ^ 

+ (e^ + 2e^ + 2e r 5 )m 2 K ml + (2e r 35 - 2e r w )m 2 K p K ■ p« + (2e r 35 - e r u ) 

mlPK ■ P* ~ &e r 39 (p K ■ pir) 2 }, (29) 



as well as 



(TT+TT-\O (27 ^^\K ) ct = ~8i^m 2 + 16i^-(4 - 2<£ + d r b - 34), (30) 



(7r+7r -|O(8,l).(l/2)|^0 )ct = m 2 + 4 -^| m 2 + 

(2eJ + 4el + 2e r n + 4e r 15 + 3e r 2 - 4e r 35 - 2e r 6 ), (31) 
for K — > 7T7T, = = m, and 



(6dl + 2^ + 2d r 20 - 8d r 2A - U\ + 2d\ + 12d r 6 ~ 9d r 7 ), (32) 



(7r+7r -| O (M),(V2)|^0 )ct = ^ (< _ m 2 } + 

(-2eJ + 6e^ + - 4^ 3 + 4e^ 5 - 4e^ - 2e^ - 4e^ + 8eJ 9 ) (33) 

for X 7r7r,m K(1 _ (oop) = 2™ 7r(1 _ Zoop) . 

From expressions (28) and (29) one can obtain the leading order LEC's 
ai, and 0:2, as well as the linear combinations e 3g ,e"[ + e 3 ,e"[ Q — e 35 + |c?g, 
and 2e\ Q — e r xl + 6dg. From Eqs. (30) and (31) for UK1 one can then obtain 
e r xl + 2e r lb — 3<ig. Making use of all the input thus obtained into Eqs. (32,33) 
for UK2 yields e r 13 — |<ig. These 14 linear combinations (namely d\, d 2 , d 7 , d\ — 

rF d r — rl r d r p r p r — p r p r A- p r p r p r — p r 4- ^rl r 9p r — p r 4- (\rl r p r 4- 
u 5 ,u 4 u 20 , u 24 , e 2 , c 1 e 5 , -t- e 3 , e 39 , e 10 e 35 t 2 «g, ze 10 L ll ^ uti 6' e ll ^ 

2ei 5 — 3dg, e^g — |dg) are sufficient to reconstruct the physical K — » 7T7t, A/ = 1/2 
amplitudes for operators such as <3i^ 2 ,<3 2 ^ 2 etc.: 
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[(-dS + d\ + 9dg - 4d?)m^ + 2(6d[ 



24 



2d r 20 



8^24 



-2dr 



6dg 



(3 



[( e ; - 2e r 13 + e\ h )m\ + {-2e\ + 2e[ + e r u + Ae[ 3 
-2e r 3 -4e r 35 + 8e r 39 )ml}. 



4e r 2 



(35) 



The logarithmic and Gasser-Leutwyler counterterm contributions to the am- 
plitudes presented in this subsection are given in Appendix D. 

4.3 (8 L ,l fl ), A/ = l/2 

The case of pure (8,1) operators, i.e. Q3,4,5,6 is simpler than the previous case 
of mixed AI = 1/2 operators; note also that phenomenologically, pure (8,l)'s 
are the most important. This is clearly a special case of the previous one for 
which (27 l, lfl) contributions are irrelevant. For the physical K — > irir reaction 
at 0(p 4 ), Eq. (35), eight new linear combinations are needed: e\~-,e\ — 65,63 + 

p r p r — p r 0p r — p r p r p r 4- 0p r p r 

e 5' e 35 6 10 ,ze 35 e n , e 39 , e n -|- ze 15 , e 13 . 

The terms quadratic in quark mass for K — > 0, Eq. (26), yield e 2 and e^-e^. 
A similar fit to K + — > 7r + , Eq. (29), then leads to eg + eg, eg 5 — e\ Q , 2e r 3i — e\ x and 
e^g. Using this for K — > 7T7t at UK1, Eq. (31), yields + 2e^ 5 and if — ■> 7T7t at 
UK2, Eq. (33), may be fitted to give e r 13 . While determining these coefficients 
is expected to be quite demanding, it is useful to note that several of them are 
obtained via more than one measurement. Note, in particular, that the term 
linear in quark mass, 02, originating from operator mixing occurs in K — ► 0, in 
K — > 7r and also in K — > 7T7t at UK1 where the operator injects energy. 

The logarithmic and Gasser-Leutwyler counterterm contributions to the am- 
plitudes presented in this subsection are a subset of those given in Appendix 
D. 



5 Conclusion 

This paper presents all of the counterterm and finite logarithm contributions to 
K° — ► |0), K° — ► K°, and K — > 7r with momentum insertion, and K — > 7T7t (at 
two values of unphysical kinematics) to 0(p 4 ) in ChPT for the (27^,1^) and 
(8l, 1_r) operators. It demonstrates that these quantities are sufficient to fully 
determine K — > tttt to 0(p 4 ) at the physical kinematics. It should be emphasized 
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that this calculation was done in full ChPT, and that these arguments do not 
necessarily apply to the quenched theory. In fact it is quite likely that some of 
the K — > 7T7T matrix elements suffer from large corrections due to the quenched 
approximation; this possibilty has recently been raised in Ref. p3[ | for the case 
of Qq. Indeed, we have done a fit to the quenched RBC data for Q^ 2 

3/2 

and Q s using the next-to-lcading order ChPT prediction of Cirigliano and 
Golowich H and have found a poor fit (x 2 /d.o.f. « 2). Thus, the data tends 
to disfavor a large coefficient for the chiral log term that is predicted by full 
ChPT. A simple quadratic fit with the coefficient of the log term set to 
yielded a much better fit (x 2 /d.o./. « 0.1). These arguments suggest that an 
unquenched lattice calculation is probably necessary in order to correctly extract 
the 0(p i ) counterterms from the lattice. It is clearly important to see whether 
this extraction procedure, especially including K — > irir at the two unphysical 
kinematics, can be extended to the case of 0(p 4 ) quenched ChPT. 

In closing we briefly want to remind the reader that two other interesting 
methods have been proposed recently for lattice extraction of K — ► ttti 

amplitudes. We believe it is important to use all the methods in order to obtain 
reliable information on this important process. 
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APPENDIX A 



Appendixes B-D contain the finite logarithm and Gasser-Leutwyler countert- 
erm 0(p 4 ) contributions to the amplitudes presented in this paper. They were 
calculated using the FeynCalc package |28| written for the Mathematica 
[ p9| system. These expressions involve the regularized Veltman-Passarino basis 
integrals A and B BQ]: 



A(m 2 ) = 



16tt 2 / 2 



2, m ,• 1 
m in — =- = hm — 



M 4 ^ 



d d £ 



(2ir) d t 2 -m 2 



2im 2 \ 



(Al) 



B(q 2 , m?, m 2 ,) = lim — 

d->4 i 



,4-d 



d d £ 



(2Tr) d ({£ + q) 2 - m()(£ 2 - m 2 2 ) 



2i\ 



13 



1 



= / dx——- [1 + ln(-x(l - x)q 2 + xm{ + (1 - x)mi) 



L "(47T) 

-ln M 2 ], 



(A2) 



where 



A = 



16tt 2 



I l I -I(ln4 7 r- 7B + l) 



(A3) 



APPENDIX B 

At 1-loop order, the pseudoscalar decay constants and masses are renormalizcd 
such that J-k.k = / (l + ^f*) and m^ (1 . lMJl) - to 2 , x (l + The 
corrections are given by 



^ = -2A(ml) - A(m 2 K ) + j^(2m 2 K + to 2 )L 4 + j^w 2 , 



^ = -^(m 2 ) - - \A{m*) + ^(2m% + to 2 )L 4 

+-^m 2 K L 5 . (B2) 



Am 2 1 1fi 

— ^ = ^ (m 2)_ A(m 2 ) + [(-L 4 + 2L 6 )2m 2 K (B3) 



-(-Li - L 5 + 2L 6 + 2L s )mi] 



Am 2 K _ 2 4/ _ 2x 16 
to 



K 



-A(m 2 ) + [(-2L 4 - L 5 + 4L 6 + 2L s )m 2 K 



(B4) 



+ (-L 4 + 2L 6 )to 2 ]. 

" :(i-; op) — "^(l-ioop) 



For degenerate quark masses at 1-loop order, m 2 K(1 _ loo , = mL_ w = to 2 ( 1 + ^ 



Am 2 9 1fim 2 

= -A(to 2 ) + ^— (-3L 4 - i 5 + 6L 6 + 2L 8 ), (B5) 
to 3 / z 

,2 



A / „ „/ 2n 8TO 2 

= -3A(to 2 ) + — 



3A(m 2 ) + — (3L 4 + L 5 ). (B6) 
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APPENDIX C 



The logarithmic corrections for the quantities relevant for the determination of 
the (27,1), AI = 3/2 K — > irn amplitudes in this paper are given by 



(K°\0%i\\K°) log 



8a 



27 



P 



16tt 2 / : 



1 13 

: m 4 K - 8m 2 K A(m 2 K ) + ( -ml - —m 2 K ) A(m 2 ) 



+ [--m 2 K - ^mlj A{ml)-2 f 



(CI) 



(7r+ | (27,l),(3/2)|^ +) 



log 



2p K ■ p 7T B(q 2 ,m 2 K ,ml) 



/ 



-A(m 2 ) - 7A(m 2 K ) 
(C2) 



(7r+ ^-| O (274),(3/2)|^0 );og 

for if — > Tnr,mK — m v — m, and 



a-" 27 2 

r 



—3m 2 ( , m 2 

5 In — + 1 



16tt 2 / 2 



3Af Am 5 
i. j 

/ 



(C3) 



(n + n-\O^>W\K ) log 



o- a 27 2 

-3i-^-m K 



K 



12 167r2/2 V 1141n-f +311,5- 1481n2 
-16cot- 1 2 + 46V^-^ 



+ 



3m 2 K 



Am 2 K - Ami 



for K — > 7T7T, ™_K-(i-( 00 p) = 2m 7I .( 1 _; oop ). The logarithmic corrections to the 
physical A/ ~ 3/2 K nn amplitude (included for completeness) are given by 



(C4) 



(n+n-\O^^\K ) log = -4i 



. «27 



L 12 



+ (rn A K - 3mlm 2 K + 2m^j B(m 2 K ,ml, ml) 
-"■i (5 + 3) A(m 2 ) + - 4< + 4< ) .!(,„* ) 



m 
4 m 2 



45 



+ 1277^ A ( m l) 



(C5) 
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where the imaginary part of expression (C5) is given by 



I m ( l (TT+n-\O (27 ' 1 ^\K )) = 



2a 27 1 



m 



(m 2 K -ml)(m 2 K -2ml). 

(C6) 



APPENDIX D 



The logarithmic corrections for the quantities relevant for the determination of 
the [(8,1)+(27,1)], AI = 1/2 K —> tttt amplitudes (as well as the pure (8,1) 
amplitudes, neglecting the (27,1) expressions) are given by 



(0\O^\K°) log 



(m 2 K - ml) 



3A(m 2 K ) 



\Aiml) 



A/. 



K 



\ (ml - 4m 2 K )A(m 2 ) - m 2 K A(m 2 K ) + \m\A{ml) 



(Dl) 



(0\O^\K°) log = ^[(ml-Am 2 K )A(ml)+4m 2 K A(m 2 K )-mlA(ml)], (D2) 



(7r +| (27,i),(i/2)|^+ )log 



27 



2 \2 



+ 7m% - §m 2 K m\ 



-m% - 36(px • p^) 2 - 6(3m 2 K - 5ml)p K ■ p^j B(q 2 ,m 2 K , m 2 ) 



+ 1 / 3(m 2 K + m\){m 2 K - ml) 2 



-3{m 2 K - ml) 2 + 20(p K ■ p^) 2 - 6{m 2 K + m\)p K ■ p n )B(q 2 ,m 2 Kl 



3 ( 7m% — 8m 2 K ml + mi. „ 2 2 \ , , 

f— * - + 9m 2 K -ml- lOp K ■ p v A(m 2 ) 



3m 2 K (m 2 K - ml 



3m 2 K + lOp K ■ Ptt A(m 2 K ) 



+ 7T 



3 / m K — m v 



/ / 



- m|- + m 2 - 6p fe • p n A(ml) 



16 



(7r+|0 (8,l),(l/2)| X+);og 



4ai 



_L_ ( (7m 2 K - ml){m 2 K - ml) 2 
72 



- 7m 4 K + 6m 2 K ml + m% 



+36(pk • V*? + 6(3m|- - 5ml)p K ■ p^j B(q 2 ,m 2 K , to 2 ) 
l('d{m 2 K + ml)(m 2 -ml) 2 



+ 



-3{m 2 K - ml) 2 + 20( PK ■ P7I ) 2 - 6(m 2 K + ml)p K ■ Pv ) B(q 2 , m 2 K , to 2 ) 
1 f -7m 4 K + 8m 2 K ml - mi 



+ - 



24 \ q 
1 / m 4 K + Am 2 K ml - 5m 4 
12 V ~Q 2 
ml 



- - 9m 2 K + ml- 30 PK ■ * A(m 2 ) 
+ \\m\ + 5ml + 30p K ■ Pit ) A(m 2 K ) 



K -m 2 K + ml- 6p K ■ Pk ) A(m%) 



Q 



Af K A.U 
f f 



PK - Pit 



2 m K 



1 / (m 2 K - ml) 



+ 



4CH2 

~¥ " I 12 \ 

1 (3(m 2 K -ml) 2 



2 \2 



-m K -m v + 6p K ■ Pn B{q , m K , m v ) 



- i{m 2 K + ml) + lOp K ■ p„ B(q 2 ,m 2 K , ml) 



_l(«pt> +7 ) 



A{ml)-\(^-/^ + 5)A{ml) 



+ - 



3 / m 2 K — m: 



--3 UK) 



2, A/ fc Af n 



f f 



(D4) 



(n+„-\OW»<W^K ) log = -8i^m 2 



for K — > 7T7T, = = m, 



-3m z 



16ir 2 f 2 



r 1 171 

51n-2- + 1 



3A/ Am 2 



/ m 2 
(D5) 



( 7r + 7r -| O (8.i),(i/2)| Jf >Jo9 



for K — > 7T7T, = m v = to, 



o- a l 2 

, - a 2 2 

r 



1 2 1 

--TO 



6"" 16tt 2 / 2 
1 



-2rrT 



16tt 2 / 2 V 9 
3A/ Am 2 32m 2 



772 

50 In ^5- -37 
M 

101, to 2 47 
In — s 

2 9 



3A/ Am 2 



TO 



(D6) 



/ 



,/' 2 



(2Li + 2L 2 + L 3 + 2i 4 + 2L 6 + L 8 ) 



(7r+ ,-| O ( 2 74),(l/2)|^0 );og = _3j^ m ^ 



-TO 



K 



24 16tt 2 / 2 V 



(456 In - 20(5 + In 16) 
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-71n5 + 232cotT 1 2 



+ 3i( Aw -- A ^ 



A/g 2A/ ff 
/ / 



(D7) 



for if -> 7T7T, m K(1 _ loop ) = 2m 7T(1 _ loop ), and 



( 7r + 7r -| O (8.1),(l/2)|j if >Iofl 



rn 



K 



72 16tt 2 / 2 
+184cot- 1 2 + 80)-^-^ 



518 In - 2091n5 + 7001n2 



-,- i Ami- — Am 2 
+12z^|m^(4L 4 - L 5 + 8L 6 + 4L 8 ), 

for if — > 7T7T, m K (Y-ioo V ) — 2m 7r ( 1 _; oop ). The logarithmic corrections to the 
physical Ai = 1/2 if — > tttt amplitude are given by 



(D8) 



(^-|0(2V),(l/2)|^O )iog = _ 4 . 



Q27 



1 B(m„,m K ,m v ) 



+ m K (7^2 - \™hc + B(m%,m 2 K ,ml) 
+ {^ 2m K + 3m 2 K ml - m^j B{m 2 K , m\,m\) 

+m 2 (m 2 K - ml^j B(m 2 K , to 2 , to 2 ) + (^f - y to 2 - + ^m 2 ^ A(to 2 ) 



c 4 

■5m 



11 



- y4 + 10m 2 ) ^(m 2 ,) 



/«4 _ 

v 2to 2 



f I ^ - 3m 2 K + -ml ) A(ml ) 



(D9) 



( 7r + 7r -| O (8.i),(i/2)| Jf >Iofl = 4 , 



2m^ - 3to^to 2 +to* )B(to^,to 2 ,to 2 ) 



1 / m^- 



-to 2 774 - to 2 B(m£-, mf l ,mf l )--[-f 



(D10) 

A(m 2 ) 
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5™^ — 3m K A(m K ) 



1 /m 4 

2 I ~ 2 

64:1 

-j E -a 2 (m 2 K - ml)[2m 2 K L 4 + (-4L 4 - £5 + 8L 6 + 4L 8 )m^] 



The imaginary parts of expressions (D9) and (D10) are given by 



Im(i(7r + Tr-\0<- 27 ^^\K Q )) = 



2a 



27 



fKfl 167T/2 



'l-4^(m%-ml)(2ml-ml), 



K 



(Dll) 



/m(i<7r+7r-|0( 8 ' 1 )'( 1 /2)|A' )) 



2ai 



16tt/ 2 



j l-A^(m 2 K -ml)(2ml-ml). 



K 



(D12) 
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